دانلود گزارشهای کارآموزی و پروژه دانشجویی

دانلود پروژه انتقال حرارت به سیالات

فهرست مطالب

فصل اول

۱-۱-سیال فوق بحرانی ۵
۱-۲-کاربردهای سیالات فوق بحرانی ۷
۱-۳-شمای کلی انتقال حرارت ۹
۱-۳-۲-انتقال حرارت در فشارهای فوق بحرانی ۱۲
۱-۳-۳-اخلال در انتقال حرارت ۱۵
شکل (۱ـ۱۳): رفتار دمای دیواره در شارهای حرارتی مختلف ۱۶
مروری بر مطالعات گذشته ۲۱
۲-۱- مقالات بازبینی ۲۲
۲-۲- مطالعات تجربی و آزمایشگاهی ۲۴
جدول ۲-۲ : تحقیقات آزمایشگاهی منتخب ۲۵
۲-۳ – تحلیل های عددی ۳۴
جدول (۲-۳) : مطالعات عددی منتخب ۳۴
۲-۴-روشهای پیش بینی
جدول (۲-۴) : روابط منتخب ضریب انتقال حرارت ۳۹
۲-۵ اخلال انتقال حرارت ۴۵
۲-۶ – اثر شتاب حرارتی ۴۹
معادلات حاکم ۵۲
۳-۱- معادلات لحظه ای حاکم : ۵۴
۳-۲- فرضیات ساده کننده : ۵۵
۳-۲- معادلات متوسط زمانی حاکم بر جریان ۵۷
۳-۴ – شرایط مرزی ۶۰
۳ـ۵ـ مدل یک بعدی ۶۱
مدلسازی و حل عددی ۶۴
۴ـ۱ـ مدل آشفتگی ۶۵
معادلات پیوستگی: ۷۷
ارزیابی مدل و بررسی نتایج ۸۳
۵ـ۱ـ پایداری حل عددی ۸۴
شکل (۵ـ۱): تأثیر تعداد مش بر نتایج حل عددی ۸۵
۵ـ۳ـ تأثیر مدل آشفتگی ۹۳
شکل (۵ـ۶): تأثیر مدل آشفتگی بر ضرایب انتقال حرارت ۹۵
۵ـ۶ـ اثر شار حرارتی ۹۹
نتیجه گیری و پیشنهادات : ۱۱۳
۶ـ۲ـپیشنهادات : ۱۱۶

فصل اول

مقدمه
انتقال حرارت به سیالات با خواص متغیر موضوعی است که از بیش از نیم قرن پیش مورد توجه محققان قرار گرفته است.
خواص ترمودینامیکی و انتقالی در سیالات معمولا تابعی از دما و فشار سیال است. این خواص در دماها و فشارهای معمولی تقریبا ثابت است. یکی از پیچیده ترین وکلی ترین سیال با خواص متغیر و تابع شدید دما و فشار سیال فوق بحرانی می باشد. این سیال بدلیل تغییر بسیار زیاد خواص آن بخصوص در نقطه بحرانی بسیار مورد توجه است و همواره به عنوان یک سیال خواص متغیر کامل مورد استفاده قرار می گیرد.
در اینجا نیز با توجه به ویژگی های این سیال که در ادامه شرح داده خواهد شد و همچنین به عنوان پیچیدهترین نوع سیال خواص متغیر که میتوان انواع دیگر از سیالات با خواص متغیر را حالت خاصی از این سیال دانست از این سیال به عنوان سیال پایه وخواص متغیر استفاده میشود.

۱-۱-سیال فوق بحرانی

وقتی صحبت از سیال فوق بحرانی میشود منظور سیال در فشار بالای نقطه بحرانی و دمای نزدیک نقطه بحرانی یا نقطه شبه بحرانی  Tpc میباشد.(شکل ۱-۱ )
نقاط شبه بحرانی به نقاطی اطلاق میشود که ظرفیت کرمایی ویژه در فشار ثابت ماکزیمم است.

شکل (۱ـ۱): نمودار درجه حرارت ـ حجم برای آب خالص

در واقع در هر فشار فوق بحرانی یک نقطه شبه بحرانی( دمای شبه بحرانی) وجود دارد که در آن تغییرات خواص سیال حداکثر است.( ظرفیت گرمای ویژه ماکزیمم است).
شکل (۱ـ۲): نمودار فشار ـ درجه حرارت برای آب خالص (دیاگرام فاز)

همانطور که در شکل (۱-۲) دیده میشود ناحیه فوق بحرانی به دو قسمت ناحیه شبه مایع  و ناحیه   شبه بخار  تقسیم میشود . در فشار ثابت زمانی که دما بزرگتر از دمای شبه بحرانی است ناحیه شبه بخار و در زمانی که کوچکتر از دمای شبه بحرانی است ناحیه شبه مایع نامیده میشود .
دلیل این اسم گذاری آن است که در واقع در فشارهای فوق بحرانی سیال را نه میتوان مایع فرض کرد ونه بخار وتنها هالهای است که فقط میتوان به آن سیال گفت . واین تقسیمبندی فقط جهت تطابق با حالت فشارها و دماهای عادی (زیر نقطه بحرانی) است وگرنه در فشارهای فوق بحرانی تغییر فاز وجود ندارد و فقط خواص سیال من جمله چگالی در قبل و بعد از نقطه شبه بحرانی تغییر میکند . همچنین میتوان اینطور عنوان کرد که بدلیل اینکه در ناحیه دمای بزرگتر از دمای شبه بحرانی چگالی کوچکتر از ناحیه دمای کوچکتر از دمای شبه بحرانی است ، ناحیه چگالی کوچکتر را شبه بخار ودیگری را شبه مایع مینامند .

۱-۲-کاربردهای سیالات فوق بحرانی
در دهه های اخیر استفاده از سیال فوق بحرانی در صنعت رو به فزونی است . برای افزایش بازده نیروگاهها در سالهای اخیر استفاده از آب فوق بحرانی SCW  ، مورد توجه قرار گرفته است . سیالات فوق بحرانی بعنوان مبرد (خنک کننده ) برای ماشینهای الکتریکی وهمچنین بعنوان مبرد برای راکتورهای هستهای مورد استفاده قرار میگیرند . در فرایندهای شیمیایی بسیار زیادی مانند ، تغییر فرم ذره  ، استخراج  و کارخانههای کف (شوینده) از سیالات فوق بحرانی مانند CO2 وهیدروکربنها استفاده میکنند .
در یکی دیگر از کاربردهای سیالات فوق بحرانی از اکسیداسیون آب فوق بحرانی ، scwo  ، استفاده میشود . مواد آلی مسموم بهمراه اکسیژن داخل آب فوق بحرانی مخلوط میشوند که نتیجه محصولات احتراق بی ضرر میباشد . در این روش آب خاصیت جالب توجهی از خود نشان میدهد. آب که در حالت معمول حلال مواد معدنی ونمکها است وقابلیت حل کردن مواد آلی را ندارد ، در حالت فوق بحرانی تغییر کرده و بر عکس میشود یعنی مواد آلی را به خوبی در خود حل میکند در حالی که دیگر حلال خوبی برای نمکهای معدنی نمیباشد . محصولات معمول از این فرایند شامل CO2 ، آب و نمکها یا اسیدهای غیرآلی میشوند .
مزایای روش scwo عبارتند از :
۱٫ تخریب سریع وکامل کربنهای آلی ، حتی با وجود دیاکسید
۲٫ عدم تولید NOx ، SOx ، دیاکسید ودوده
۳٫ محصولات تولیدی آب ، دیاکسید کربن و گاز نیتروژن است .
۴٫ بازیافت انرژی قابل حصول است .
شکل (۱-۳): کثافات فاضلاب
(سمت چپ:پیش از اعمال،سمت راست:پس از اعمال)

روش scwo در موارد ذیل قابل اعمال است :
آبهای آلوده آلی شامل ؛ آبهای پسرفتی کارخانههای شیمیایی ، آبهای پسرفتی صنایع غذایی
کثافات آلی شامل ؛ کثافات فاضلاب شهرداری ، کثافات کارخانههای شیمیایی
در بیشتر این کاربردها ، انتقال حرارت به فرایند سیال در فشارهای فوق بحرانی ودر دماهای زیر نقطه بحرانی و فوق بحرانی وجود دارد .

۱-۳-شمای کلی انتقال حرارت
۱-۳-۱-خواص فیزیکی حرارتی
انتقال حرارت در فشار فوق بحرانی بشدت نشات گرفته از خواص فیزیکی حرارتی است که بشدت تغییر میکنند (بخصوص در نزدیکی خط شبه بحرانی).
شکل (۱-۴ ) نماینگر رابطه بین ظرفیت گرمایی ویژه Cp وفشار ودما میباشد .
شکل (۱ـ۴): ظرفیت گرمایی آب

همانطور که دیده میشود در هر فشار یک ظرفیت گرمایی ویژه ماکزیمم محلی وجود دارد. در ناحیه فشار زیر بحرانی  حداکثر مقدار ظرفیت گرمایی ویژه بر روی خط اشباع قرار میگیرد . در نقطه بحرانی ( C۳۷۴ =T و MPa1/22=P ) ظرفیت گرمایی ویژه بیشترین مقدار خود را دارد .
در فشارهای فوق بحرانی ، مکان هندسی نقاطی که مقدار ماکزیمم ظرفیتهای گرمایی ویژه را به هم وصل میکند خط شبه بحرانی ، PCL  ، نامیده میشود که در شکل(۱-۵) نشان داده شده است .
در فشار  MPa25 دمای نقطه شبه بحرانی  C384 است . ظرفیت گرمایی ویژه در نقطه بحرانی kJ/kg K5600 (شکل(۱- ۶)) است که ۱۰۰۰ برابر بزرگتر از این مقدار در دمای معمولی است .
اشکال (۱- ۷ ) تا (۱-۱۰)  نمایانگر تغییرات خواص فیزیکی حرارتی در برابر دما در فشارهای مختلف میباشد .
شکل(۱ـ۶): گرمای ویژه در PCL                      شکل (۱ـ۵): خط شبه بحرانی PCL در یک دیاگرام P-T
شکل (۱ـ۸): هدایت گرمایی آب فوق بحرانی                      شکل (۱ـ۷): چگالی آب فوق بحرانی
شکل (۱ـ۱۰): عدد پرانتل آب فوق بحرانی                      شکل (۱ـ۹): ویسکوزیته آب فوق بحرانی
در نزدیک خط شبه بحرانی (PCL ) ، چگالی  بشدت کاهش مییابد . در این ناحیه ضریب انبساط حرارتی دارای پیک بزرگی است که رفتاری شبیه ظرفیت گرمایی ویژه از خود نشان میدهد .ضریب هدایت حرارتی  با افزایش دما کاهش مییابد ، هر چند در نزدیکی نقطه شبه بحرانی یک ماکزیمم محلی وجود دارد . ضریب هدایت حرارتی پس از دمای شبه بحرانی با شیب تندی کاهش مییابد . ویسکوزیته دینامیکی  نیز رفتار مشابهی از خود نشان میدهد . در نقطه شبه بحرانی برای عدد پرانتل  بدلیل افزایش  ظرفیت گرمایی ویژه با شیب تند ، پیک بزرگی اتفاق میافتد .

۱-۳-۲-انتقال حرارت در فشارهای فوق بحرانی
همانطور که در بخش قبل گفته شد ، در نزدیکی خط شبه بحرانی تغییرات بسیاری در خواص فیزیکی حرارتی اتفاق میافتد . این قضیه میتواند باعث تغییرات بسیار شدید در ضریب انتقال حرارت گردد . با محاسبه معادله دیتوس-بولتر
(۱ـ۲)

شکل (۱ـ۱۱): ضریب انتقال حرارت با توجه به معادله دیتوس ـ بولتر

برای جریان آب ، در حالت آشفتگی در یک لوله مدور ، واستفاده از دمای حجمی  جهت محاسبه خواص در آن حالت ، ضریب انتقال حرارت بدست میآید . در شکل (۱-۱۱) ضریب انتقال حرارت بر حسب دمای حجمی در شار جرمی   Mg/m2s1/1 ، فشار  MPa25 ،شار حرارتی   MW/m28/0 و قطر لوله  mm4 نشان داده شده است.
همانطور که دیده میشود در نقطه شبه بحرانی (  ˚ C384 =T ) معادله دیتوس-بولتر ضریب انتقال حرارت را برابر kW/m2K40 میدهد که تقریباً بیش از دو برابر مقدار آن در دماهای پایین  (برای مثال ˚C300) و پنج برابر مقدار آن در دماهای بالاتر (برای مثال ˚C500) است . این قضیه به وضوح نشان میدهد که بدلیل تغییرات در خواص فیزیکی حرارتی ، ضریب انتقال حرارت در نزدیکی خط شبه بحرانی به شدت تغییر میکند . هر چه فشار به فشار نقطه بحرانی نزدیکتر باشد ، پیک ضریب انتقال حرارت بلندتر میشود .
همچنین در این اشکال دیده میشود که ضریب انتقال حرارت که توسط معادله دیتوس-بولتر بدست میآید از مقدار واقعی آن ، بخصوص نزدیک خط شبه بحرانی ، انحراف دارد . در شارهای حرارتی کوچک ، ضریب انتقال حرارت از مقداری که بوسیله  معادله دیتوس-بولتر تخمین زده میشود بزرگتر است . این پدیده افزایش انتقال حرارت  نامیده میشود . در شارهای حرارتی بزرگ، ضریب انتقال حرارت از مقداری که توسط معادله دیتوس-بولتر بدست میآید کوچکتر است . شکل(۱-۱۲) نسبت ضریب انتقال حرارت  را به مقدار محاسبه شده توسط معادله دیتوس-بولتر  را نشان می دهد . همانطور که دیده می شود تحت شرایطی خاص این نسبت بسیار کوچک می شود .
شکل (۱ـ۱۲): نسبت ضریب انتقال حرارت   به مقدار محاسبه شده با معادله (۱ـ۱)

۱-۳-۳-اخلال در انتقال حرارت
انتقال حرارت در یک لوله را در نظر بگیرید . در شکل(۱-۱۳) دمای دیواره  لوله بر حسب دمای حجمی سیال رسم شده است . منحنی ها به ترتیب نمایانگر هر دو حالت شار حرارتی کوچک و شار حرارتی بزرگ هستند . در شار حرارتی کوچک دمای دیواره رفتار ملایم و یکنواختی را نشان میدهد و با افزایش دمای حجمی افزایش مییابد . در این حالت اختلاف بین دمای دیواره و دمای حجمی کوچک باقی میماند .
در یک شار حرارتی بزرگ نیز رفتار مشابهی ، به جز در زمانی که دمای حجمی به مقدار شبه بحرانی نزدیک میشود ، از دمای دیواره مشاهده میشود .در این حالت افزایشی با شیب تند در دمای دیواره اتفاق میافتد . وقتی دمای حجمی از دمای شبه بحرانی میگذرد ، دمای دیواره دوباره کاهش مییابد . این افزایش بسیار زیاد در دمای دیواره مربوط به پدیده “ اخلال در انتقال حرارت” میشود . هنوز در مقالات مختلف تعریف واحدی برای شروع اخلال ر انتقال حرارت وجود ندارد. تعدادی از مقالات بحران جوشش را به عنوان دلیل این پدیده میدانند . این در حالی است که کاهش ضریب انتقال حرارت ، یا افزایش دمای دیواره در مقایسه با رفتار بحران جوشش  ، که در آن  دمای دیواره  با شیب تندی افزایش مییابد ، رفتاری ملایمتر ویکنواختتراز خود نشان میدهد.
شکل (۱ـ۱۳): رفتار دمای دیواره در شارهای حرارتی مختلف

اخلال در شارهای حرارتی بزرگ نسبت به شار جرمی،  G/ًq، ایجاد میشود. این اخلال ناشی از دو مقوله جدا از هم میباشد. یکی بدلیل تغییر ساختار آشفتگی مرتبط با گرانش (g) در نزدیکی دیواره و دیگری بدلیل اثر دمپینگ آشفتگی  ناشی از شتاب حرارتی میباشد. اولی را اثر شناوری  و دومی را اثر شتاب حرارتی  مینامند.
الف ـ  اثر شناوری
در جریانهای سیالات با خواص ثابت، وجه غالب در انتقال گرما، جابجایی اجباری میباشد. جریانهای سیالات با خواص متغیر (شرایط فوق بحرانی) اتقال حرارت بصورت ترکیبی از جابجایی آزاد و اجباری خواهدبود. جابجای آزاد، در واقع ناشی از تأثیر نیروهای شناوری است. در مورد لوله های قائم، شناوری موجب میشود که بین ضرایب انتقال حرارت بدست آمده برای  جریان بالارو و پایینرو اختلاف اساسی ایجاد شود. انتقال حرارت در جریان پایینرو افزایش مییابد در حالیکه درجرانهای بالارو انتقال حرارت دچار اخلال و زوال میگردد. در لولههای افقی در شرایطی که قطر لوله نسبتاً بزرگ، شار جرمی  کم و شار حرارتی زیاد باشد اثرات شناوری نمود میکند. در این لولهها اثرات شناوری باعث تغییرات محیطی در ضرایب انتقال حرارت میشود، بطوریکه مقادیر ضریب انتقال حرارت در قسمتهای فوقانی لوله کمتر از مقادیر مربوطه در قسمتهای تحتانی لوله میشود. اختلال ایجاد شده در انتقال حرارت در قسمت فوقانی لولههای افقی را میتوان به لایهبندی و طبقهبندی  در جریان نسبت داد. در واقع در لولهها افقی در حالتی که شار جرمی کم، شار حرارتی زیاد و دمای دیواره از دمای شبه بحرانی بیشتر شده باشد، سیال در نزدیکی دیواره به حالت شبه گاز و در قسمت مرکزی لوله به حالت شبه مایع در می آید .در این حالت دمای لایه های نزدیک دیواره به دمای شبه بحرانی رسیده و باعث تغییرات بسیار زیاد در خواص سیال میشود. یکی از این تغییرات  کاهش چشمگیر چگالی میباشد. لذا به همین دلیل بخش شبه مایع سیال که چگالی بیشتر دارد جای خود را با بخش شبه گاز در دیواره پایینی لوله که چگالی کمتری دارد عوض میکند. بدین ترتب در نزدیک دیواره پایینی لوله، سیال شبه مایع و در نزدیک دیواره فوقانی لوله سیال شبه گاز قرار میگیرد. از آنجا که ضریب هدایت حرارتی بخش شبه مایع بزرگتر است، پس فرایند انتقال حرارت در قسمت تحتانی لوله نسبت به قسمت فوقانی لوله بهتر انجام میگیرد.

(( برای دانلود کلیک کنید ))